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We investigate here spherically symmetric gravitational collapse in a spacetime with an arbitrary 
number of dimensions and with a general type I matter field, which is a broad class that includes 
most of the physically reasonable matter forms. We show that given the initial data for matter in 
terms of the initial density and pressure profiles at an initial surface t — U from which the collapse 
evolves, there exist rest of the initial data functions and classes of solutions of Einstein equations 
which we construct here, such that the spacetime evolution goes to a final state which is either a 
black hole or a naked singularity, depending on the nature of initial data and evolutions chosen, and 
subject to validity of the weak energy condition. The results are discussed and analyzed in the light 
of the cosmic censorship hypothesis in black hole physics. The formalism here combines the earlier 
results on gravitational collapse in four dimensions in a unified treatment. Also the earlier work is 
generalized to higher dimensional spacetimes to allow a study of the effect of number of dimensions 
on the possible final outcome of the collapse in terms of either a black hole or naked singularity. 
No restriction is adopted on the number of dimensions, and other limiting assumptions such as 
self-similarity of spacetime are avoided, in order to keep the treatment general. Our methodology 
allows to consider to an extent the genericity and stability aspects related to the occurrence of naked 
singularities in gravitational collapse. 
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I. INTRODUCTION 
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A considerable amount of work has continued in recent years on the validity or otherwise of the cosmic censorship 
conjecture (CCC) in black hole physics. The reason for this interest is that CCC is fundamental to many aspects of 
the basic theory and applications of black holes, and the astrophysical implications resulting from the phenomena of 
^ ■ continual gravitational collapse of a massive star which has exhausted its nuclear fuel. As of today, no theoretical 
I | proof, or even any satisfactory mathematical formulation of CCC is available, where as an intriguing finding that has 
Vh . emerged from recent investigations on gravitational collapse scenarios in general relativity is that the final end state 
of such a collapse could be either a black hole (BH), or a naked singularity (NS). The theoretical and observational 
J> . properties of these objects could be quite different from each other and so it is of much interest to get an insight into 
how each of these phases come about as end states of a dynamically developing collapse governed by gravitational 
dynamics. 

When a sufficiently massive star starts collapsing gravitationally on exhausting its nuclear fuel, it would not settle 
to a stable configuration such as a neutron star. What happens in such a case is an endless gravitational collapse 
ensues, where the sole governing force is gravity. According to general relativity, the outcome of such a process would 
be necessarily a spacetime singularity. If an event horizon of gravity forms well in advance before the singularity 
forms that gives rise to a black hole as the final state of collapse. In the case otherwise, when the horizon formation is 
delayed during the collapse, these extreme density and curvature regions may fail to be covered by the horizon, and a 
visible naked singularity may develop (see e.g. [l[ - 0] for some recent reviews). The possible physical consequences 
of the later scenario have drawn some attention recently [l[ . Most of the gravitational collapse models studied so far 
are spherical, and the matter cloud collapses under reasonable physical conditions such as an energy condition, and 
regularity conditions on the initial data from which the collapse develops. 

A note-worthy suggestion that has emerged towards a possible theoretical formulation of CCC is that, any naked 
singularities resulting from matter models which may also develop singularities in special relativity, should not be 
regarded as physical 0] - [l2j ■ Clearly, it will require a serious effort to cast this into a mathematical statement and 
a possible proof for CCC. Also, it may not be easy to discard completely all the matter fields such as dust, perfect 
fluids, and matter with various other reasonable equations of state, which have been studied and used extensively in 
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relativistic astrophysics for a long time from the perspective of understanding gravitational collapse processes and 
final states. 

Another possibility that indeed appears worth exploring is we may actually be living in a higher dimensional (HD) 
spacetime. The recent developments in string theory and other related field theories strongly indicate that gravity 
is possibly a higher dimensional interaction, which reduces to the general relativistic description at lower energies. 
Hence, there is a possibility that while CCC may possibly fail in the four-dimensional manifold of general relativity, 
it may well be restored in higher dimensions due to the extra physical effects arising from our transition itself to a 
higher-dimensional spacetime continuum. Such considerations would inspire a study of gravitational collapse in higher 
dimensional spacetimes. From such a perspective, many works have reported results in recent years on spherically 
symmetric collapse in HD. The recent revival of interest in this problem is partly motivated by the Randall-Sundrum 
brane-world scenario [l3[ and different 'sectors' of the general problem have been studied by restric ting to certain 
subcases such as those in self-similar spacetimes, or spacetimes with a specific number of dimensions |l4j - (22j. 

It follows that a general investigation of gravitational collapse in a spacetime with an arbitrary number of higher 
dimensions will be of considerable interest. From the perspective of CCC, the effects the number of dimensions may 
have on the final outcome of gravitational collapse in terms of BH/NS phases will be of much interest. 

From such a perspective, we investigate here the issue of BH/NS endstate formations in a higher dimensional 
gravitational collapse in some detail, in order to bring out in a transparent manner the effect of dimensions on the 
final fate of evolution of a matter cloud which collapses from a given regular matter initial data. A spherically 
symmetric collapse is considered here in N > 3 dimensions and the matter content is chosen to be of Type I, which 
obeys the weak energy condition. The collapsing matter field we consider here is a general and broad class, which 
includes most of the physically reasonable matter fields such as dust, perfect fluids, massless scalar fields and such 
others and restrictions of any special form are not imposed on the form of matter. 

We consider spherically symmetric spacetimes here, however, when compared to some of the earlier studies of 
collapse in higher dimensions mentioned above, we deal here with the general case. That is, there is no further 
restriction imposed on the number of dimensions, the spacetime is not assumed to obey various special conditions 
such as self-similarity or assuming the existence of various Killing fields, which are somewhat restrictive assumptions 
for collapse models. The results here also generalize earlier results [H|-[2(|, [27|, related to gravitational collapse 
final states, and provide a unified treatment. 

In order to investigate collapse final states, given the initial data for matter in terms of the initial density and 
pressure profiles at an initial surface t = ti from which the collapse develops, we construct classes of solutions to 
Einstein equations, such that the spacetime evolution goes to a final state which is either a black hole or naked 
singularity, depending on the nature of rest of the free initial data functions and possible evolutions, subject to 
validity of the weak energy condition. The methodology used here allows us to consider in some detail the genericity 
and stability aspects related to the occurrence of naked singularities in gravitational collapse. 

In the next Section II, we describe the basic equations and also the various regularity conditions for gravitational 
collapse. Section III then considers spherical collapsing clouds and their basic dynamics. The apparent horizon 
and structure of trapped surface formation is discussed in Section IV, which has a direct bearing on the nature of 
singularity in terms of being either visible or covered. The issues related to equation of state and that of validity 
of energy conditions in our consideration are discussed in Section V, and we also discuss briefly collapse in lower 
spacetime dimensions in Section IV. The exterior spacetime and matching conditions are given in Section VII and the 
final Section summarizes some concluding remarks. 

II. EINSTEIN EQUATIONS, REGULARITY AND ENERGY CONDITIONS 

Let us consider a general N-dimensional spherically symmetric metric of the form, 

ds 2 = -g ab (x",x 1 )dx a dx b + R 2 (x", x l )dn 2 N _ 2 (1) 



where a, b run from to 1, and 
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is the metric on (N — 2) sphere with 0* being the spherical coordinates. From this metric, we get the elements of 
Einstein tensor as [Hj [lj| , 

N — 2 

N G ab = [g ab {2RUR + (N - 3)Q} - 2RR >ab ] (3) 



3 



N G 22 = -i [i? 2 5ft + (N- 3){2RDR + (N - 4)Q}] (4) 



N G u {i>2) 

where we have, 
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G22 (5) 



Q= l + i?- a i? a ; Di? = g ab R. ab (6) 

and, 

5K = g ab M ab (7) 

Here 5R a fc is the Ricci tensor evaluated by the two-metric g a b, and 5ft is the scalar curvature evaluated by the same. Also, 
throughout this paper we use the ususal convention, that is, a comma(,) in subscript denotes partial differentiation 
while a semicolon (;) denotes covariant differentiation. 

Let us now describe the spacetime geometry within the spherically symmetric collapsing cloud by the comoving 
coordinates (t,r,9 l ), which are specified as below. We take the matter field to be of Type I, which is a broad class 
including most of the physically reasonable matter forms, including dust, perfect fluids, massless scalar fields and such 
others. This is specified by the requirement that the energy momentum tensor for the matter admits one timelike and 
(N — 1) spacelike eigenvectors [2^|. We now choose our co-ordinates (t,r,9 l ) to be along these eigenvectors, which 
makes the coordinate system to be comoving, that is, the co-ordinate system moves with the matter. We can use the 
freedom of coordinate transformations of the form t' = f(t, r) and r' = g(t, r) to make the gt r term in metric (fTJ) and 
the radial velocity of the matter to vanish. In that case the general metric in the comoving coordinates (t, r, 9 1 ) must 
have three arbitrary functions of t and r and this can be written in the form (29j . 

ds 2 = _ e Mt,r) dt 2 + e 2i,(t,r) dr 2 + R 2 ^ (8) 

In this comoving frame the energy-momentum tensor for any matter field which is Type I is given in a diagonal form, 

7? = -p(t,r); T r r =p r {t,r); T$=p e (t,r) (9) 

The quantities p, p r and pg are respectively the energy density, and radial and tangential pressures, ascribed to the 
matter field. The matter cloud has a compact support with < r < rb, where Tb denotes the boundary of the cloud, 
outside which it is to be suitably matched through suitable junction conditions with another spacetime geometry. 

We take the matter field to satisfy the weak energy condition, that is, the energy density measured by any local 
timelike observer is non-negative. This ensures the physical reasonability for the collapsing matter fields we are 
considering. Another energy condition frequently used is the dominant energy condition, which demands that for any 
timelike observer the local energy flow is non-spacelike. We note that these two are frequently regarded as the main 
and important energy conditions which are physically reasonable. For these energy conditions to be satisfied, we 
must have for any timelike vector V 1 , 

T lk V l V k > (10) 
and T ik Vk non-spacelike. For the energy-momentum tensor @, these amount respectively to the conditions, 

P>0; p + p r >0; p+p e >0 (11) 

\pA<p ;\pe\<P- (12) 

Now with the above metric (8), the following quantities can be evaluated, 

R a R. a = -R 2 e- 2v + R 2 e- 2 ^ (13) 

OR = -e~ 2v [R + R(-v + J;)] + e- 2 ^'[R" + R'(v' - V')] (14) 
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and, 

R ab = (r-»a- < -^-) sy b + (sr - ip'R 1 - sy b + Q ab (u) 

where we have, 

Q ab = (# - v'R - i>&) M + 615°) (16) 

Using these quantities the Einstein equations Gik = Tik take the form, (in the units 8irG = c = 1), 

_ (N-2)F< (N-2)F 
P ~ 2R N ~ 2 R> ' Pr ~ 2R N ~2R [ ' 

v ,_ (N-2)(po-p r ) R p' r 

p+Pr R P+ Pr 

- 2R! + R% + R?L = (19) 



where we have defined, 

G(t, r) = e- 2 ^{R'f- H(t, r) = e- 2v {R) 2 (21) 

The function F = F(t, r) is known as the Misner Sharp mass, which gives the total mass in a shell of comoving 
radius r, at an epoch t. The energy condition p > imply F > and F' > 0. Since the area radius vanishes at the 
center of the cloud, from equation (|17[) it is evident that in order to preserve the regularity of density and pressures 
at any non-singular epoch t, we must have -F(t,0) = 0, that is the mass function should vanish at the center of the 
cloud. 

As seen from equation (fTT)) . there is a density singularity in the spacetime at R = 0, and at R = 0. However, the 
later ones are due to shell-crossings [30j , which basically indicates the breakdown of the coordinate system we have 
used. These are not generally regarded as genuine singularities, which can be possibly removed from the spacetime to 
extend the manifold through the same [3l|. Hence we shall consider here only the shell- focusing singularity at R = 0, 
which is a genuine physical singularity where all matter shells collapse to a zero physical radius. We shall discuss this 
in some more detail in the next section. 

We note that, in general, for a general matter field with non-vanishing pressures as we consider here, there are a 
variety of dynamical time evolutions possible from the given matter density and pressure profiles as prescribed on an 
initial surface (which we call here matter initial data), from which the collapse evolves. In particular, even if the cloud 
commences gravitational collapse at the initial surface t = U, there can be classes of solutions of Einstein equations 
where the evolution is such that a bounce is possible at a later stage for the cloud. We consider here only continually 
collapsing class of models, because our interest is in the physical situation which corresponds to the case when the 
mass of the star is so high that on exhausting its nuclear fuel, it must undergo a continual gravitational collapse, 
completing the same in a finite time. Thus the continual collapse condition is included as a part of our work here. In 
this case, trapped surfaces develop and a spacetime singularity necessarily forms as collapse end state and we need to 
find the conditions when the final singularity is necessarily covered within an event horizon (as hypothesized by the 
cosmic censorship), or when it will be naked with strong gravity regions being visible to faraway observers. In the 
case of a bounce or dispersal, no singularity of course need to form in the spacetime, a situation with which we do 
not concern ourselves here. 

We can use the scaling freedom available for the radial co-ordinate r to write R = r at the initial epoch t = tj. It 
is interesting to note that, if we had wished to scale the radial co-ordinate at the initial epoch as R(ti,r) — r 13 , (3 
being any constant, then the only possible allowed value for (3 is unity. Because in the case otherwise, either R would 
blow up at the center (which is not allowed as the Einstein's equations would require the metric functions to be at 
least C 2 ), or R' would go to zero at the center causing a shell-crossing singularity (which we would like to avoid, by 
construction), which violates the regularity of the initial data. 
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We now introduce a function v(t, r) as denned by, 

v(t,r) = R/r (22) 

We then have R(t, r) — rv(t, r), and 

«(t(,r) = l; «(t,(r),r)=0; £>< (23) 

The time t — t s (r), that is v = 0, corresponds to the shell-focusing singularity at R — 0, which is a genuine spacetime 
singularity where all the matter shells collapse to a vanishing physical radius. The condition v < here corresponds 
to a continual collapse of the cloud. The description of shell-focusing singularity at R = in terms of the function 
v(t, r) has several advantages. The physical radius goes to the zero value at the shell-focusing singularity, but we also 
have R = at the regular center of the cloud at r = 0. This is to be distinguished from the genuine singularity at 
the collapse end state by the fact, for example, that the density and other physical quantities including the curvature 
scalars all remain finite at the regular center r = 0, even though R = holds there. This is achieved, as we point out 
below, by a suitable behavior of the mass function, which should go to a vanishing value sufficiently fast in the limit 
of approach to the regular center where (even though R goes to zero) the density must remain finite. On the other 
hand, when we use the function v(t, r), we note that at t = i$ we have v = 1 on the entire initial surface, and then as 
the collapse evolves, the function v continuously decreases to become zero only at the singularity t s {r), that is, v = 
uniquely corresponds to the genuine spacetime singularity at R = 0. 

^From the point of view of dynamic evolution of the initial data prescribed at the initial epoch t — U, there are five 
arbitrary functions of the comoving shell-radius r [26J , as given by, 

v(ti,r) = p (r), ip(t h r) = ip (r), R(t h r)=r, (24) 
p(U,r) = p (r), p r (U,r) =p ro (r), pe(U,r) = pe 

We note that not all the initial data above are mutually independent, because from equation (fl~8|) we get, 

Mr)= n ( "- 2) ^ (25, 

Jo V r{po+p ro ) P0+Pr J 

Thus, apart from the matter initial data describing the initial density and pressure profiles, the rest of the initial data 
which is free is ipo(r), which essentially describes the velocities of the collapsing matter shells as we shall discuss later. 

To ensure regularity of the initial data, the initial pressures must be taken to have physically reasonable behavior 
at the center. Considering that the total force at the center of the collapsing cloud should be zero, we have the 
gradients of initial pressures vanishing at the center. Also, regularity of the initial data requires that we must have, 
JV o (0) — f>0o(O) = 0' that is, at the center the difference between the radial and tangential pressure vanishes. The 
metric functions have to be C 2 differentiable everywhere as per the requirements of the Einstein equations, and as 
seen from the equation for z/, the above condition is implied by the requirement that v' does not blow at the regular 
center. This means that the matter should behave like a perfect fluid at the center of the cloud with the net force 
vanishing there. We note that these regularity conditions do not exclude the purely tangential pressure (p r = 0) or 
purely radial pressure collapse models (pg = 0) which we shall refer to later, because in those cases the above implies 
that pe — > or p r — > respectively, close to the center, where the matter then closely approximates dust. We note 
that these regularity conditions give us a sufficient condition for the regularity of the metric function vq(t) at any 
non-singular initial epoch. It follows from equation (f2"5|) that at the center of the cloud both vq and v' go to zero. 
Hence vq{t) has the form, 

vo{r) =r 2 g(r) (26) 

where, g(r) is an arbitrary function which is at least C 1 for r = 0, and it is at least a C 2 function for r > 0, as Einstein 
equations demand the metric functions to be at least C 2 everywhere. Another regularity condition frequently used in 
collapse considerations is there are no trapped surfaces at the initial surface from which the collapse begins. 

We thus see that there are five total field equations with seven unknowns, p, p r , pe, ip, v, R, and F, giving us the 
freedom of choice of two free functions. Selection of these free functions, subject to the weak energy condition and 
the given regular initial data for collapse at the initial surface, determines the matter distribution and metric of the 
space-time, and thus leads to a particular time evolution of the initial matter and velocity distributions. As we shall 
show, it turns out that given the matter initial profiles in terms of pa,p ro and pg , there exist rest of the initial data 
at t — ti, and classes of solutions, which we find by means of explicit construction, which give either a black hole or 
a naked singularity as the end state of collapse. The outcome depends on the nature of rest of the initial functions, 
and the classes of dynamical evolutions as allowed by the Einstein equations. 
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An important point to be noted here is, in the description above we have made no mention so far on the equation 
of state that the matter must obey. Typically, these are of the form, p r — p r {p) and pg — pe{p). If these are specified, 
then there is no freedom left, and we have seven equations for seven variables. If we are to incorporate this right 
away, the only way to proceed to find the collapse end state would be to assume a specific equation of state that the 
matter must satisfy, and then to examine the collapse problem and the nature of the final singularity as resulting from 
the dynamical evolution as governed by the Einstein equations. There have been many collapse studies in past using 
such an approach, e.g. for dust equation of state, perfect fluids etc. The limitation of such an approach, however, has 
been that there is very little existing knowledge on what a realistic equation of state should be that the matter has 
to satisfy at the extreme high densities that a continual collapse realizes in its advanced stages. For example, even 
for neutron star densities which are relatively low as compared to those of continual collapse, there is a great deal of 
uncertainty on the equation of state for such neutron matter. As a result, the neutron star mass limits are uncertain 
to that extent. Thus, specific or special assumptions used on the equation of state may turn out to be physically 
unrealistic or restrictive and untenable in the final stages of collapse. In fact, diametrically opposite views exist on 
the possible equation of state in very late stages of collapse. For example, while there are many arguments suggesting 
that pressures must play important role in the later stages of collapse, the opposite view is that in such late stages 
the matter must necessarily be dustlike (see e.g. (32|, (33|) . 

Under the situation, the path we take here is, we do not assume any specific or particular equation of state 
presently, and carry out our further considerations in a general way in terms of the allowed initial matter profiles, 
and the allowed dynamical evolutions of the Einstein equations, to determine the black hole and naked singularity 
end states for collapse. We then discuss subsequently, in Section V, the role that the equation of state will play 
towards further fine tuning the BH/NS outcomes as collapse end states. The advantage such an approach has is, 
first we write all the collapse equations in generality, and then only different subcases can be examined depending on 
the corresponding equation of state under consideration. As we shall show, various important subcases such as dust, 
perfect fluids and others can be included as special cases of the treatment given here. 

In this paper our basic strategy is as described below. We actually do not choose any explicit form of F(t,r), it is 
allowed to be completely general class, subject to regularity conditions. Till Section V, we do not make any special 
choices of the two free functions, that we refer to, but deal with and construct certain general classes, as determined 
by the differentiability conditions on the concerned functions (e.g. as specified by the equation (45) onwards) so that 
basically the singularity curve (which is defined as the time taken for a shell labelled V to reach the singularuty) is 
expandable upto at least first order, i.e. the singularity curve is regular enough so as to have a well-defined tangent, 
for these classes of evolutions. Then we show that these classes contains both BH and NS final states, as decided by 
the sign of the tangent of singularity curve. 

In other words, given the matter initial data, we construct rest of the functions so that in the continual collapse 
(R < 0), the final singularity curve is regular (expandable) as above with a well-defined tangent. With that we then 
show, when the final singularity is visible, or covered in a black hole. The key point is, once we have the matter initial 
data, we show the existence of classes of rest of the initial data and evolutions which are solutions to the Einstein 
equations, by explicit construction as specified here, so that the final state is either a BH or NS depending on the 
choice made of the rest of the initial data and evolutions. This is subject to energy conditions and other regularity 
conditions. We are of course not dealing with all possible classes of evolutions from the given matter initial data, 
which is not our purpsoe, but we just show the existence of the classes that lead the collapse to above final states. 

For a given matter initial data, once we have chosen the rest of the initial data, and the evolutions as solutions to 
the Einstein equations, that clearly fixes the equation of state for the matter. That can be quite 'exotic' at times, 
depending on the choices made. However, as shown explicitly in Section V, the construction given here does include 
well-known equations of state, such as dust, perfect fluids etc. 

III. GRAVITATIONAL COLLAPSE OF MATTER CLOUDS 

The method we follow is outlined as below. In the case of a black hole developing as collapse end state, the spacetime 
singularity is necessarily hidden behind the event horizon of gravity, whereas in the case of a naked singularity 
developing there are families of future directed non-spacelike trajectories which terminate in the past at the singularity, 
which can in principle communicate information to faraway observers in the spacetime. The existence of such families 
confirms the NS formation, as opposed to a BH collapse end state. We study the singularity curve produced by the 
collapsing matter, and it is shown that the tangent to the same at the central singularity at r = is related to the 
radially outgoing null geodesies from the singularity, if there are any. By determining the nature of the singularity 
curve and its relation to the initial data and the classes of collapse evolutions, we are able to deduce whether the 
trapped surface formation in collapse takes place before or after the singularity. It is this causal structure of the 
trapped region forming during collapse that determines the possible emergence or otherwise of non-spacelike curves 
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from the singularity. This settles the final outcome of collapse in terms of either a BH or NS. 

Given the matter initial profiles in terms of the functions Po{r) , p ro (r) , p@ (r) at the initial epoch t = i, from which 
the collapse commences, our purpose now is to construct and examine possible evolutions (classes of solutions to 
Einstein equations) of such a matter cloud to investigate its final states. 

While constructing the classes of solutions which give the collapse evolutions, given the matter initial data at t = tj, 
we preserve as much generality as possible. Hence we allow the mass function F(t, r) for the collapsing cloud to have 
a general form as given by, 



F(t,r) = r"- l M(r,v) 



(27) 



where M. > is at least a C 1 function of V for r = 0, and at least a C 2 function for r > 0. 

It is to be noted that F must have this general form follows from the regularity and finiteness of the density profile 
at the initial epoch t — U, and at all other later regular epochs before the cloud collapses to the final singularity at 
R = 0. This requires, from the Einstein equation (17), that F must behave as r^ N ~^ close to the regular center. 
Hence we note that since M. is a general (at least C 2 ) function, the equation equation (j2"T|) is not really any ansatz 
or a special choice, but quite a generic class of the mass profiles for the collapsing cloud, consistent with and allowed 
by the regularity conditions. We thus make no special choice of F but allow it to be a general function as given by 
the above equation. Then equation (I17| gives, 



p(r,v) 



N -2 



(N-l)M + r[M, r + M, v v r \ 



and 



p r (r,v) = - 



(AT -2) 



„JV-2 



The regular density distribution at the initial epoch is given by, 



Po(r) = ^77^ [(N - l)M(r, 1) + rM(r, l), r 



(28) 



(29) 



(30) 



It is evident that, in general, as v — ► 0, p — ► oo and p r — > oo. That is, both the density and radial pressure blow up 
at the shell-focusing singularity. 

It is seen that given any regular initial density and pressure profiles for the matter cloud from which the collapse 
develops, there always exist energy profiles or velocity functions for the collapsing matter shells, and classes of 
dynamical evolutions as determined by the Einstein equations, so that the collapse end state would be either a naked 
singularity or a black hole, depending on nature of the allowed choice. Thus, given the matter initial data at the 
initial surface t = U, these evolutions take the collapse to end either as a BH or NS depending on the choice of the 
class, subject to regularity and energy conditions. 

To see this, we need to construct classes of solutions to Einstein equations to this effect. Let us define a suitably 
diffcrentiable function A(r, v) as follows, 



u'(r,v) = A{r,v), v R' 



(31) 



That is, A(r,v). v = v' /R', and since at t = ti we have R = r which gives [A(r, v) tV ] x = v' (r). Our main interest 
here is in studying the shell-focusing singularity at R = which is the physical singularity where all the matter shells 
collapse to zero radius. Therefore we assume that there are no shell-crossing singularities in the spacetime where 
R' = 0, and that the function A(r, v) is well-defined. From equation (|26)l . we generalize and choose the form of v(t, r) 
as the class given by, 



v(t,r) = r 2 gi(r,v) 



(32) 



where g\ (r, v) is a suitably diffcrentiable function and g\ (r, 1) = g(r). It then follows that A(r, v) has the form, 

A(r,v) = rg 2 (r,v) (33) 

From the regularity conditions and that the total force at the center of the collapsing cloud should be zero at any 
non-singular epoch, it is evident that both gi(r, v) and g%(r,v) shoud be well defined at r — and v ^ 0. In fact in 
some of the models discussed in section V, these functions are well behaved even at the singularity. In general, it is 
possible that these functions may in fact blow up at the singularity in some cases. But, as we prove later, even if 
these functions do blow up, the singularity curve can still be well defined and expandable. 
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Some comments are in order on our assumption that R' > 0, that is, we have considered here the situation with 
no shell-crossing singularities. This is because, it is generally believed (sec e.g. [30] 31]) that such singularities can 
be possibly removed from the spacetime as they are typically gravitationally weak, and because spacetime extensions 
have been constructed through the same in certain cases [341 ]. In contrast, in several physically reasonable collapse 
models including dust and perfect fluids R = turns out to be a gravitationally strong curvature singularity. Under 
the situation, we are interested only in examining the nature of the shell- focusing singularities at R = 0, which are 
genuine curvature singularities arising at the termination of collapse, where the physical radii for all collapsing shells 
vanish, and the spacetime necessarily terminates without extension. 

Specifically, R' > implies that we must have v + rv' > 0. Since v is necessarily positive through out the collapse, it 
follow that this will be satisfied always whenever v 1 is greater or equal to zero. Even when it is negative, the condition 
that the magnitude of rv' should be less then that of v is sufficient to ensure that there will be no shell-crosses. Later 
in this section we derive an expression for the quantity v', in terms of the initial data and the other free evolutions 
as allowed by the Einstein equations. Hence it follows that we can specifically state the condition for avoidance of 
shell-crossings in terms of the behavior of these functions. 

Coming to the dynamical collapse evolutions, using equation (l3T|) in equation (I19|) we get, as a class of solutions of 
Einstein's equations, 

G{r,v) = b(r)e 2rA{r ' v) (34) 

Here b(r) is another arbitrary function of the shell radius r. By the regularity condition on the function v at the 
center of the cloud we get the form of b(r) as, 

b(r) = 1 + r 2 b (r) (35) 
where bo(r) is the energy distribution function for the shells. Using equation (|3ip in equation (| 18[) . we get, 

Rv' 

(N - 2) Pe = RA tV (p + Pr ) + (N- 2) Pr + (36) 

In general, both the density and radial pressure blow up at the singularity, so the above equation implies that the 
tangential pressure would also typically blow up at the singularity. Now using equations l|27p. (f3Tj) and ((34)) in 
equation ([20"]) , we get, 



R^R = -e< r ^ + r 2 b Q )R( N - 3 )e 2rA ^ v 1 - i?^" 3 ) + A N ~ X )M (37) 

The negative sign in the RHS of the above equations corresponds to a collapse scenario where we have R < 0. Defining 
a function h(r,v) as, 



the equation (f20|) becomes 



h(r, v) = 1 = 2g 2 (r, v) +0(r 2 ) (38) 



v^v = -y^+^M^-^o + & 2v (v<- N - 3 )h + M) (39) 
Integrating the above equation with respect to v, we get, 

t(v,r)= ■ (40) 



« ^Je( rA +^v( N -Vbo + e 2v {v^ N -^h + M) 



Note that the variable r is treated as a constant in the above equation. The above equation gives the time taken for 
a shell labeled r to reach a particular epoch v from the initial epoch v = 1. Expanding t(v, r) around the center of 
the cloud, we get, 

t(v,r) =t(v,0)+rX(v) + 0(r 2 ) (41) 

where the function X(v) is given as, 

_ _i ^,"(,(»»^M») + A<,(,)) 

2j„ (vf-<>)b K + v<. N -Vli ll +Mo(v))t 
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wherein we have defined, 

b OQ = b Q (0); M o (v)=M(0,v); h = h(0,v) 

&i = 6 (G); M 1 (v)=M, r (0,v); hi = h r (0,v) (43) 

Hence we see that the time taken for a shell labeled r to reach the spacetime singularity at R = (which is the 
singularity curve) is given as, 

/ n / V 2 CIV r i t\ 

t s (r) = / , — (44) 

JO J e (rA+u) v (N-3) bo + e 2u (y(N-3) h + M ) 

As we want to consider here continual collapse, we focus only on those classes of solutions where t s (r) is finite and 
sufficiently regular. This means that the cloud collapses in a finite amount of time. In the physical situation of a 
continual collapse of a massive matter cloud in a finite amount of time the function t s (r) has to be of course finite. As 
for regularity, to check the existence conditions for a well defined, continuous and C 2 singularity curve, let us define 
a function Q(r, v) as, 

N — 3 

Q(r, v) = . (45) 

J e (rA+u) v (N-3) bo + e 1v ( v (N-3)h + M) 

Also we consider the following functions, 

4>i(r) = / Q{r, v)^ r dv ; <j> 2 (r) = / Q(r,v), rr dv (46) 



Jo Jo 
Let A be the rectangular area in (r, v) plane defined by the lines, 

r = 0;r = e ;v = 0;v = 1 (47) 
Now if the following conditions are satisfied [35[ , 
f . Q(r, v) is a continuous function of r and v in A, 

2. Q(r, ?j). r and Q{r, v)^ rr are continuous functions of r and t> in A, 

3. The integrals 01(7") and fair) converge uniformly in A. 
then we can write, 

Mr) = A. [ts{r)] . Mr) = JP_ [ts{r)] (48) 

and this implies that the singularity curve t 3 (r) would be a well-defined C 2 function near the center. 

We would like to emphasize here that regularity of different functions that makes Q{r, v), namely is, A, and M, is 
sufficient but by no means necessary condition for the existance of a C 2 singularity curve. As we can easily see, even 
if these functions blow up at r = 0, v = 0, the function Q may still be well defined. For example, if we suppose these 
functions tend to infinity at the singularity the function Q would tend to zero. For any regular collapse, the function 
t s (r) has to be finite and well-defined. 

Several well-studied collapse models such as dust collapse and others satisfy these or stronger conditions, and so 
the singularity curve is well-defined and expandable. Some examples of such singularity curves are as follows. In case 
of an N-dimensional dust collapse, the expression of the singularity curve would be [361 ] . 

t. r ) dust = / —==_==== (49) 

where A4(r) and bo(r) are well defined C 2 functions of the comoving coordinate r and well defined at r = 0. ft can 
then be easily seen that the singularity curve is differentiable at the center, with, 

dt s (r) _ I r 1 v^jMi+v^-^b^dv 
dr ~ 2 Jo (Mo+v( N - 3 )b 00 )i 
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where we have defined, 

6 00 = &o(0), M =M(0) 
h =b'(0), M 1 =M'(0) (51) 

Another such example is an N-dimensional Einstein Cluster, which describes a non-steady spherically symmetrical 
system of non-colliding particles moving in such a way that relative to a suitably moving frame of co-ordinates, their 
motion is purely transversal. In this case, the singularity curve is given by [37j . 



v 2 \ v z -\ — V- civ 



ts(r) Ec = / (52) 

Here L(r) is a function of the radial coordinate r only. In the Newtonian limit, this function corresponds to the angular 
momentum per unit mass of the system. Therefore the function L(r) is called as the specific angular momentum, and 
has the form, 

L(r) = r 2 l{r) (53) 

Again, since A4(r) and bo(r) and L(r) are well-defined C 2 functions of the coordinate r, we see that the above 
singularity curve is well defined and differentiable at r = 

To generalize this and to give another explicit example, it can be shown that given any initial data of the form 
(|2"5)) . there always exist classes of dynamical evolutions, which give rise to a well defined and differentiable singularity 
curve. Let us, by the freedom of choice of free functions, choose the evolution functions M.(r, v) and A(r,v) in the 
following way, 

M(r, v) = m(r) — p r (r)v N ~ 1 , A(r, v\ v = va(x\ Xi x = rv (54) 

From Einstein equations it is clear that for the above class of evolutions the radial pressure remains static. However, 
the tangential pressure blows up along with the density at the singularity v — and is given by, 

2 Pe (r,v) = p r +p' r (R/R') + v {rv) tR [p(r,v)+p r ] (55) 

One can now easily check that the above class of evolutions admits a well-defined and differentiable singularity curve 
because both the functions M.(r, v) and A(r, v) are well-defined and C 2 at r — and v = 0. 

Once we have a singularity curve which is at least C 2 , we can Taylor expand the function near the center as, 

t s (r) =t So +rX(0) + O(r 2 ) (56) 

where t so is the time when the central singularity at R = 0, r = develops, and is given as, 

From the above equation it is clear that for t so to be defined, 

v^ N - 3) b 00 + v^ N ^h +M (v) >0 (58) 



In other words, a continual collapse in finite time ensures that the above condition holds. Also, from equation (J39J) 
and (141 [) we get for small values of r, along constant v surfaces, 

v^v' = sj (v( N - 3 )b 00 + v( N -Vh + M (v))X(v) + 0{r) (59) 

It is now clear that the value of X(0) depends on the functions bo, M. and h, which in turn depend on the initial data 
at t — ti, the dynamical variable v, and the evolution function A(r, v). Thus, a given set of initial matter distributions 
and the dynamical profiles including the energy distribution of shells completely determine the tangent at the center 
to the singularity curve. 
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IV. APPARENT HORIZON AND THE NATURE OF THE SINGULARITY 



It is now possible to examine, given the matter initial data at the initial surface t = tt, how the final fate of collapse 
is determined in terms of either a black hole or a naked singularity. If there are families of future directed non-spacelike 
trajectories reaching faraway observers in spacetime, which terminate in the past at the singularity, then we have a 
naked singularity forming as the collapse final state and in the case otherwise when no such families exist and event 
horizon forms sufficiently earlier than the singularity to cover it, we have a black hole. This is decided by the causal 
behavior of the trapped surfaces developing in the spacetime during the collapse evolution, and the apparent horizon, 
which is the boundary of trapped surface region in the spacetime. 

In general, the equation of apparent horizon in a spherically symmetric spacetime is given as, 

u' h i:.,i:.L o (60) 

Thus we see that at the boundary of the trapped region the vector R ± is null. Substituting ([5]) in ([50]) we get, 

R' 2 e" 2 ^ - R 2 e~ 2v = (61) 
Using equation (|20[) we can now write the equation of apparent horizon as, 

F 



R N - 



1 (62) 



which gives the boundary of the trapped surface region of the spacetime. If the neighborhood of the center gets 
trapped prior to the epoch of singularity, then it is covered and a black hole results, otherwise it could be naked when 
non-spacelike future directed trajectories escape from it. 

Thus the important point is to determine if there are any future-directed non-spacelike paths emerging from the 
singularity. To investigate this, and to examine the nature of the central singularity at R = 0, r = 0, let us consider 
the equation for outgoing radial null geodesies which is given by, 



dt 
dr 



(63) 



We want to examine if there would be any families of future directed null geodesies coming out of the singularity, thus 
causing a naked singularity phase as collapse endstate. The singularity occurs at v(t s (r),r) = 0, i.e. R(t s (r),r) = 0. 
Therefore, if there are any future directed null geodesies terminating in the past at the singularity, we must have 
R — > as t — » t s along these curves. Now writing equation ([63]) in terms of variables (u — r a ,R) , we have, 



— = -r-^-^R' 
du a 



(64) 




In order to get the expression of the tangent to null geodesies emerging in the (R, u) plane, we choose a particular 
value of a such that the geodesic equation is expressed only in terms of known limits. For example if X(0) ^ 0, and 
the functions M. and h are well defined for < r < rt and < v < 1 we choose a — t^4- Using equation (|20p and 

considering R < 0, we then get the null geodesic equation in the form, 

_ F \ 

(65) 

If the null geodesies do terminate at the singularity in the past with a definite tangent, then at the singularity the 
tangent to the geodesies have 4^ > in the (u, R) plane, and must have a finite value. In the case of a massive 
singularity in dimensions greater than or equal to four,(ie. F(t s (r),r) > for r 0), all singularities for r > are 
covered since r n-b ~^ °o and hence ^ — > — oo. This is when both the pressures p r and pg are positive with the 
energy condition being satisfied. Therefore in such a case only the central singularity at R — 0, r = could be naked. 

Hence we need to examine the central singularity at r — 0, R = to determine if it is visible or not and to 
determine if there are any solutions existing to the outgoing null geodesies equation, which terminate in the past 
at the singularity, going to faraway observers in future. The conditions are to be determined under which this can 
happen. We also note that, since the singularity curve and the evolution functions are regular, we can calculate the 
limit of the functions H, G and F/R at r — > 0, t — > t So . From equation (|34[) . as A(r, v) is a well defined function, we 
can see that G(t So ,0) — 1. Also from equation (|39p we see that at this point H w r 2 /v N ~ 3 . Calculating this limit 
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on t = t So plane from equation (|59|) . at the point (t So ,0), we have H = 0. Hence we see from equation (|20|) that 
F/R N - 3 = in this limit. 

Let now a?o be the tangent to the outgoing null geodesies in (R, u) plane, at the central singularity, then it is given 
by, 



,. ,. R dR 
Xo = lim hm — = — — 

t— ft, r^o u du 



(66) 



t— >t.:r— >0 



To find out whether the null geodesic equation admits any solution of xq which is positive and finite at the central 
singularity, we can use the values of H, G and F/R at (t so , 0) in equation (f6"5)) . Also we use equation to get the 



value of v'v ? , on v — surface at 



(that is, on the point (t So ,0)). Thus solving equation (f6"5|) , we get, 

N-l N-l 



■y/Mo(fl)X(0) 



(67) 



and the equation of radial null geodesic emerging from the singularity is given by R = xqu in the (i?, u) plane, or in 
(t, r) coordinates it is given by 



V + 1 

t-t s (0) =x r«-i 



(68) 



It follows now that if A'(O) > 0, then xo > 0, and we get radially outgoing null geodesies coming out from the 
singularity, giving rise to a naked central singularity. However, if X(Q) < we have a black hole solution, as there will 
be no such trajectories coming out. If X(0) = then we will have to take into account the next higher order non-zero 
term in the singularity curve equation, and a similar analysis has to be carried out by choosing a different value of a. 

To show that the above is a necessary as well as sufficient condition for an outgoing radial null geodesic emerging 
form the singularity to exist, let us assume that such geodesies do exist and in the (i?, u) plane, it's equation is 
R = xqu, and xo > 0. Then at the central singularity (R = 0, u — 0), the tangent to such geodesic must be xq. Also 
this tangent must be the root of the equation, 



dR 
du 



N- 1 




n + i y u (|)^ 

at the point (R = 0, u = 0). This is possible if and only if 

'N-l 




= 



(69) 



Xq = 



v/M (0)A-(0) 



(70) 



and for the slope to be defined and positive we must have X(0) > 0. 

Now to see that X(0) > is a sufficient condition for the existence of an outgoing radial null geodesic emerging from 
the singularity, let us consider the case that the singularity curve has a positive tangent at the central singularity. 
Consider now the curve, 



t-t s (0) 



N-l 



y6Vfo(0)Af(0) 



N + l 
J- N-l 



(71) 



Along this curve t — ► t so as r — > 0. And as we have X(0) > 0, this curve is outgoing in the sense that t increases as 
we increase r along the curve. Let us now calculate the quantity {—gttdt 2 + g rr dr 2 ) along this curve in the vicinity of 
the central singularity. Using (f59|) . (f39|) and ([34]) . we have for this curve at the point (i So ,0), 



e 2v dt 2 + e~ 2rA R' 2 dr 2 = 



N + l 

N - 1 



N - 1 



-i N-l 



v/-M (0)*(0) 



~(-dr 2 +dr 2 ) = 



(72) 



That is, in the vicinity of the central singularity the curve (|71[) is null. Thus we see that given any positive value of the 
tangent to the singularity curve at the central singularity, we can always find a null and outgoing curve terminating 
in the past at the central singularity, making the singularity naked. 

We make below some remarks on the nature of the apparent horizon and its relation with the visibility or otherwise 
of the singularity. To find the equation of apparent horizon near the central singularity, let the time corresponding to 
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a shell labeled by r entering the apparent horizon, in term of the variable v, be v a h(r). Then from equation (|62[) . we 
can easily see that v a h{r) is the root of the equation, 

r 2 M{r,v)-v (N -V =0 (73) 

Now using equation (|40p , we get the the equation for apparent horizon in (t, r) plane as, 

t ah {r) = f.(r) - / (74) 

Jo J e (rA+ V ) v (N-3) bo + e 2u (y(N-3) h + M ) 

It is obvious that the necessary condition for the existence of a locally naked singularity is that the apparent horizon 
curve must be an increasing function at the central singularity, in the lowest power of r. 

We note that in the above the functions h and M. are expanded with respect to r around r = and the first-order 
terms are considered. At times, however, these are assumed to be expandable with respect to r 2 , and it is argued 
that such smooth functions would be physically more relevant. Such an assumption comes from the analyticity with 
respect to the local Minkowskian coordinates (see e.g. 38]), and it is really the freedom of definition mathematically. 
We may remark that the formalism as discussed above would work for such smooth functions also, which is a special 
case of the above discussion. 

We thus see how the initial data in terms of the free functions available determine the BH/NS phases as collapse 
end states, because X(Q) is determined by these initial and dynamical profiles as given by equation (|4"!?|) . It is clear, 
therefore, that given any regular initial density and pressure profiles for the matter cloud from which the collapse 
develops, we can always choose velocity profiles so that the end state of the collapse would be either a naked singularity 
or a black hole, and vice-versa. 

We believe numerical work on collapse models may provide further insights into these interesting dynamical phe- 
nomena, and especially when collapse is non-spherical, which remains a major open problem to be considered [4l|-[53j. 
Numerical and some analytical works have been done in recent years on spherical scalar field collapse [54l]-|62j] and 
also on some perfect fluid models [63|-[7lj]. While we have worked out explicitly here the emergence of null geodesies 
from the singularity, thus showing it to be naked (or otherwise) in an analytic manner, the numerical simulations 
generally discuss the formation or otherwise of trapped surfaces and apparent horizon, and such considerations may 
possibly break down closer to the epoch of actual singularity formation. In that case, this may not allow for actual 
detection of BH/NS end states, whereas important insights on critical phenomena and dispersal have already been 
gained through numerical methods. Probably, a detailed numerical investigation of the structure of null geodesies in 
collapse models may provide further insights here. 



V. EQUATION OF STATE AND ENERGY CONDITIONS 

As stated above, we work here with type I matter fields, which is a rather general form of matter. However, it is 
important to note that suitable care must be taken in interpreting these results. While we have shown that the initial 
data and dynamical evolutions chosen do determine the BH/NS end states for collapse, the point is, actually, all these 
dynamical variables are not explicitly determined by the initial data given at the initial epoch (note that v plays the 
role of a time coordinate here). Hence these functions are fully determined only as a result of time development of 
the system from the initial data provided we have the relation between the density and pressures, that is a given 
'equation of state'. 

In principle, it is possible to choose these functions freely (e.g. the matter and velocity profiles at the initial epoch 
and the dynamical evolutions such as F(v,r) and v(v, r)), only subject to an energy condition and regularity, which 
then fully determines the collapse evolution. One can then calculate the energy density, and the radial and tangential 
pressures for the matter. However, in that case, the resultant 'equation of state' could be quite strange in general. If 
any equation of state of the form p r = /(/?) and/or pg — g(p) is given, then it is clear from equation (|17|) and (|18|) 
that there would be a constraint on the otherwise arbitrary function M. and A, specifying the required class, if the 
solution of the constraint equation exists. It is certainly true that, presently we have practically very little idea on 
what kind of an equation of state should the matter follow, especially at very high densities and closer to the collapse 
end states, where we are already dealing with ultra-high energies and pressures. Hence if we allow for the possibility 
that we could freely choose the property of the matter fields as above, or the equation of state, then our analysis is 
certainly valid and give several useful conclusions on possible collapse end states. In such a case, it is also possible 
that the chosen equation of state will be in general such that the pressures may explicitly depend not only on the 
energy density, but also on the time coordinate. 
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All the same, it is important to point out that the analysis as given above in fact does include several well-known 
equations of state and useful classes of collapse models, also satisfying the energy conditions throughout the collapse, 
as we demonstrate below. 



A. Dust collapse 

The idealized class of dust collapse models where the pressures are taken to be vanishing has been studied extensively 
so far and has yielded many important insights on collapse evolutions. In this special case, the Einstein equations 
can be solved completely to get the N-dimensional generalization of the usual Tolman-Bondi-Lemaitre (TBL) dust 
collapse metric ([72|-[74j]) and it is given as, 

ds2 = dt2 ~ i i R L( ^ - R2 ^ r ) d ^N-2 (75) 
1 + Hoo(r) 



The equations of motion are given by, 



and 



(N - 2)F' 

2Rl»-»R> = 9 (?6) 



# = JS) + m (77) 

In the case of dust, the mass function must be F = F(r), and hence regularity condition implies that, 

F(r)=r^ N -^M(r) (78) 

The energy condition here gives < r < ?*(,, and so we must have A4(r) > and 3M + rM. r > 0. In this case, as we 
have already seen in section III, the function X(y) is given as, 



and the time taken for the central shell to reach the singularity is given by, 

iV-3 

v 2 dv 



y/Mo + v( N -Vb 



(80) 



It is now seen clearly that any given sets of density and velocity profiles at the initial epoch completely determine the 
tangent to the singularity curve at the central singularity. Also the equation (|67p becomes 

N - 1 

(81) 

It therefore follows that, given any specific density profile of the collapsing dust cloud, we can always choose a 
velocity profile so that the end state of the collapse would be either a naked singularity or a black hole depending on 
the choice made, such that energy conditions are satisfied throughout the collapse. The converse also holds, namely 
one can choose a given velocity profile for the cloud at the initial epoch, and then there are density profiles which 
will lead the collapse to either of the BH/NS final states, and these conclusions hold irrespective of the number of 
dimensions of the spacetime. Hence we see that our results unify and generalize the earlier results of dust collapse 
[7^|-[8l|. Basically the point that follows here is that, given an initial density profile for the collapsing cloud, the 
space of velocity profile functions is divided into the regions that lead the collapse cither to a black hole or naked 
singularity evolution, depending on the choice made, and the converse holds similarly. 



B. Collapse with static radial pressure 



While the dust equation of state discussed above is fairly standard and extensively used, it is widely believed that 
pressures could play an important role in gravitational collapse considerations. We discuss below a class of collapse 
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models with non-zero pressures, which is however idealized in the sense that while the tangential pressure can be 
arbitrary, the radial pressure is taken to be static. 

As we already pointed out in section III, if we consider the classes of collapse in which the radial pressure remains 
static, the constraint equation for M. has the following solution, 

M(r,v) = m{r)-p r {r)v N - 1 (82) 

In addition to this if there is an equation of state of the form pg = f(p), then that gives the constraint equation for 
the function A(r, v) in the following way. 

V(p) =Pr+ p' r (R/R') + A(r, v), v \p(r, v) + pr] (83) 

For this class of models the energy conditions are given by, 

(N - l)[m - p r v N ^\ + r[m, r - Pr ir v N ~ X - (N - l)p r v N ~ 2 v'] > (84) 

(N - l)[m - p^- 1 ] + r[m, r - p r r u Ar - 1 — (N — l)p r v N - 2 v'] + p r R 2 R' > (85) 



2 2 



,R 
R' 



> (86) 



As shown in [82J there exist classes of functions m, p r and A such that naked singularity is the end state for the 
collapse and also the above three energy conditions are satisfied. It follows that, given initial matter profiles, there 
exist classes of collapse evolutions satisfying the energy conditions as we see above, such that either of the BH/NS 
endstates can result subject to above equation of state. 

C. Isentropic perfect fluid with a linear equation of state 

Perfect fluids have been widely used in astrophysical considerations and a linear equation of state is well-studied. 
We discuss below how the formalism outlined here apply to this case to find the BH/NS configurations as a perfect 
fluid collapse end states. 

For an isentropic perfect fluid, whose pressure is a linear function of the density only, the equation of state of the 
collapsing matter is given by, 

p r (t,r) =p e (t,r) =kp(t,r) (87) 

where k € [—1, 1] is a constant. The case k — gives the dust case we discussed above and k = 1 is the stiff fluid case. 
Let us at present consider only the case of positive pressures. In that case k > and the energy conditions give, 

M,„ < (88) 

^From the above equation of state and the Einstein equations we can immediately see that the function M. is now 
the solution of the equation 

(N - l)kM + krM, r + Q(r, v)M, v = (89) 

where, 

Q(r,v) = (k + l)rv' + v (90) 
Now the above equation (f8"9"|) has a general solution of the form (83[, 

T(X,Y) = (91) 
where X(r,v,A4) and Y(r,v,A4) are the solutions of the system of equations, 

dM dr dv 



(N-l)kM kr 



(92) 
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Thus we can easily see that equation (|89|) admits classes of solutions when v' > 0. Also solving the equation for the 
central shell r — 0, with boundary conditions p — > oo as v — ► 0, we get, 

By choosing too > 0, we can make the central shell to satisfy the energy condition p(t, 0) > for all epochs. Then 
by the continuity of the density function, we can say that there exists an e-ball around the central shell for which 
v'(t,r) > and also p(t,r) > 0. But as know, at the central singularity y/vv' ~ X(Q), hence this implies that we can 
have classes of solutions which satisfy the energy conditions and also admits a naked singularity as the collapse end 
state. 

For further discussion on perfect fluid collapse, and the details of black hole and naked singularity formation, we 
refer to (63|- [7l[, and references there in. 

It is seen from the above that several well-known classes of collapse models form subcases of the consideration given 
here. Along with these well-known models, the above analysis would work for any other models with other equations 
of state, if that permit solutions to the constraint equations on M. and A. Hence it follows that the considerations 
above provide an interesting framework for the study of dynamical collapse, which is one of the most important open 
problems in gravity physics today. 

VI. THE CASE OF 2+1 DIMENSIONAL COLLAPSE 

Several studies on gravitational collapse scenarios in (2 + 1) dimensional spacetimes, have been carried out by 
various authors [Hj]- [13] • These provide interesting toy models which may provide quite important insights from 
the perspective of quantum gravity. This is because in (2 + 1) dimensional spacetimes there is no gravity outside 
matter. Also the spacetime metric is always conformally flat as the Weyl tensor vanishes identically everywhere. The 
situation in (3 + 1) and higher dimensions is far more complicated as compared to this. 

To investigate the final outcome in (2 + 1) dimensional collapse, let us consider the geometry of the trapped surfaces 
in this case in some detail. From equation (|62[) we see that the equation of apparent horizon in this case is given by, 

F(t,r) = l (94) 

It is interesting to note that the geometry of trapped surfaces here is completely determined by the mass function 
of the cloud, and is independent of the area radius of the collapsing shells. If the mass function of the collapsing 
configuration is bounded from above, say for example with F(t,r) < 1 for t £ [— oo, t 8 {r)] 1 and r e [0, r&] (where r;, 
is the boundary of the collapsing cloud) , we then see that the trapping does not occur and the complete singularity 
that forms as collapse end state is necessarily visible to an outside observer. This is strikingly different from four or 
higher dimensional cases where a massive singularity is always trapped. 

An interesting subcase of this situation is that of (2 + 1) dimensional dust collapse, when F(t,r) = F(r) with 
the mass function having no time dependence. Here we see that the initial mass of the collapsing cloud completely 
determines the final outcome in terms of BH or NS. Clouds with small enough mass always form a visible singularity, 
whereas for larger masses a trapped region is present at all epochs. However, as demanded by the regularity conditions, 
we must avoid trapped surfaces on the initial surface t = t, where the collapse commences. In that case, for the dust 
collapse case, there are no trapped surfaces developing at all at any other later epochs till the singularity formation, 
and as a result the (2 + 1) dust collapse always necessarily produces a visible naked singularity, as opposed to BH/NS 
phases obtained in usual four dimensional dust collapse which we discussed earlier. 

VII. EXTERIOR SPACETIME AND MATCHING CONDITIONS 

To complete a collapse model, we need to match the interior spacetime to a suitable exterior spacetime. As we are 
interested here in modelling collapse of astrophysical objects (such as massive stars), we have assumed the matter to 
have compact support at the initial surface, with the boundary of the cloud being at some r — r&. If one assumes the 
pressures at the boundary of the cloud to be vanishing, then it is always possible to match the interior spacetime with 
an empty Schwarzschild exterior. However, in all cases it may not be possible to make the pressures at the boundary 
of the cloud to vanish. Hence we outline here the procedure to match the interior with a general class of exterior 
metrics, which are the generalized Vaidya spacetimes [88[ and I89J. at the boundary hypersurface £ given by r = Tb- 
For the required matching we use the Israel-Darmois conditions ([9ll| and [92|). where we match the first and second 
fundamental forms, (the metric coefficients and the extrinsic curvature respectively) at the boundary of the cloud. 
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Whereas the procedures used below are standard, we shall describe the particular case treated here in some detail 
so as to give the exact picture of the overall collapse scenario emerging. We note a useful fact that since we are 
matc hing the second fundamental form Kij , there is no surface stress energy or surface tension at the boundary (see 
e.g. [93]). The metric just inside X is, 

dg 2_ = _ e 2u(t,r) dt 2 + e mt,r) dr 2 + R 2 ^ ^2 (Q5) 

which describes the geometry of the collapsing cloud. The metric in the exterior of X is given by, 

ds\ = - (l - 2M( ^' V) j dV 2 - 2dVdr v + r 2 v dn 2 (96) 

where V is the retarded (exploding) null co-ordinate and r v is the Vaidya radius. Matching the area radius at the 
boundary we get, 

R(r b ,t) = r v (V) (97) 
Then on the hypersurface X, the interior and exterior metrics are given by, 



and 



ds 2 s _ = - e 2 »^dt 2 + R 2 (t,r b ) 2 dCl 2 



<'4 + = -[i- 2M{r r v,V) +^)<n 2 >-:,i<f 



Matching the first fundamental form gives, 



dt ) y _ /-, 2M(r„V) , 2<i~ 
r v dV 



(r ) s = R(t,n) (ioo) 



Next, to match the second fundamental forms (extrinsic curvatures) for the interior and exterior metrics, we note 
that the normal to the hypersurface X, as calculated from the interior metric, is given as, 



and the non-vanishing components of the normal as derived from the generalized Vaidya spacetime are, 



(101) 



nX = 1 (102) 

'-. _ 2M(r„V) , o dr„ 
r v dV 



j _ 2M(r„,V) drv_ 
n r = ^ — jjL_ (103) 

L 2M(r„V) pdr„ 
v v dV 



Here the extrinsic curvature is defined as, 



K ab = \c n g ab (104) 



That is, the second fundamental form is the Lie derivative of the metric with respect to the normal vector n. The 
above equation is equivalent to, 

K ab = ~ [gab, c n c + g c bn c a + g ac n c b \ (105) 
Now setting \Kqq — Kg~ g ] s = on the hypersurface X we get, 

^ _ 2M(r v ,V) dr^ 

RR'e^ = r v r " =JL= (106) 

L 2M(r v ,V) 2 rfr„ 
r v dV 
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Simplifying the above equation using equation (|100[) and the Einstein equations, we get, 

F(t,r b ) = 2M(r v ,V) (107) 
Using the above equation and ()100|) we now get, 

fdV\ e^+jQ 

[~dt ~ (108) 

Finally, setting [if,^ — ^^Is = 0, where r is the proper time on S, we get, 

M(r v , V), rv = — + —f=-VH, f + Re 2 "v'e-^ (109) 

Any generalized Vaidya mass function M(v,r v ), which satisfies equation (|109|) will then give a unique exterior space- 
time with required equations of motion given by other matching conditions, (|107|) . (|108|) and |97|) . 

To see that the set of all such functions M(v, r v ) is non-empty, we have the examples of a charged Vaidya spacetime 
M = M (V) + Q(V)/r v , and the anisotropic de-Sitter spacetime M = M(r v ) as two different solutions of the equation 
(I109P (see for example [89| and |90(). This gives two unique exterior spacetimes, both of which are subclasses of the 
generalized Vaidya metric. 



VIII. CONCLUDING REMARKS 



We make here some remarks towards a conclusion and note some unresolved and open issues. 

1. Towards investigating endstates of a continual gravitational collapse, given a general type I matter field and 
given the matter initial profiles at the initial surface from which the collapse develops, we constructed here classes 
of solutions to the Einstein equations such that the collapse evolution goes to the formation of either of a black 
hole or naked singularity endstate, depending on the choice of the evolution made and choice of rest of the initial 
data functions such as velocities of the collapsing shells. This is subject to satisfying energy condition, and several 
reasonable and important equations of state are included in the framework here. It also becomes clear that in higher 
dimensions also both black holes and naked singularities can occur as collapse endstates. 

2. We like to note that what we have deduced here is the occurrence of a locally naked singularity only, as opposed 
to that of a globally naked singularity. That is, we show when the null geodesies escape from the spacetime singularity, 
going out in the future, but we do not address the question of when they go out of the boundary of the matter cloud. 
It is possible, in principle, that the singularity is only locally naked and trajectories do come out but they all fall back 
into the singularity again, without going out of the boundary of the star, thus not being globally visible. 

This issue is still not studied for the general class of models such as the classes we considered here and may be of 
interest. However, for the case of dust models this has been studied in some detail and it is shown that whenever the 
singularity is locally naked, one can always choose the classes of the mass and energy functions suitably, as one moves 
away from the center, in such a manner that the singularity becomes globally visible. The point is, while the local 
visibility of the central singularity is basically decided by the conditions near the center, the global visibility really 
depends on the overall behavior of these functions within the matter cloud, away from the center. This we are still 
free to choose. In other words, for the dust collapse models once the singularity is locally visible, there are always 
classes of functions which we can choose so as to make it globally visible. 

Another important related point here is, as such there is no scale in the problem, and the size of the collapsing 
cloud could be quite large. In such a case, even if the singularity is only locally visible, still it can be seen for a long 
enough time by the observers. Thus, in principle, a locally naked singularity is also as serious violation of the cosmic 
censorship as a globally visible singularity, and there may not be a qualitative difference in the two cases in many 
situations of physical interest. 

3. Another important and interesting issue frequently mentioned regarding occurrence of naked singularities in 
gravitational collapse is their genericity and stability. It is argued that if these are not generic or stable, these need 
not be taken seriously. This is interesting because in general relativity there is no well-defined notion or criteria 
available for stability or genericity, which one can then apply and test for a given model to ascertain these. On the 
other hand, a consideration of this issue is important all the same in that, depending on the collapse situation under 
discussion one would like to formulate these notions in some way to examine if the naked singularities developing as 
collapse endstates are 'generic' or 'stable' in some suitable sense. This would typically involve taking into account the 
topology and metric of the concerned function spaces which define the given collapse scenario. Without discussing 
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this here in further detail, we note that since we have formulated here the collapse in a general manner, for general 
physically reasonable matter fields, this classes constructed here may provide a good arena to explore and test these 
important issues for the cosmic censorship hypothesis. 

4. A related issue of course is that of non-spherical collapse. Our considerations here are restricted to spherical 
collapse, and the question remains open regarding the final end state of a non-spherical collapse. So one can ask if 
the conclusions available in spherically symmetric collapse remain the same and stable under possible non-spherical 
perturbations. 

Even though some non-spherical collapse models have been discussed and investigated so far such as the Szekeres 
quasi-spherical collapse or some cylindrical collapse models [4l|- (53|, these are somewhat restricted in nature. A 
strong theorem about the formation of trapped surfaces in cylindrically symmetric spacetimes is given in 52 1 . 
Recently, a numerical construction of naked singular solutions with the cylinderical symmetry was done in 53 1. 



It is not clear as yet if this issue can be approached really in an analytic manner, and possibly detailed numerical 
simulations of the collapsing stars could be the answer. 

5. A question that is frequently asked in connection to the occurrence of naked singularities as collapse endstates 
is that how to understand this phenomena physically. A naked singularity signifies the escape of light and particle 
trajectories from the ultra-dense spacetime regions. However, the gravity must become so strong in these regions. In 
such a case, how can any thing escape at all from such a region is the question. Thus, while a black hole which is 
a region from which not even light would escape, may appear to be the only physically reasonable outcome in such 
situations, formation of a naked singularity in collapse may appear to be counter-intuitive. 

The point that comes out from considerations such as ours is that the naked singularities are more an artifact of 
general relativity, rather than that of a purely Newtonian physics. Even though the matter density grows higher and 
higher without bound and blows up closer to a spacetime singularity, which would denote the growth of attractive 
forces of gravity, there are other important factors which are purely general relativistic effects which can delay the 
formation of trapped surfaces governing the trapping of light. 

An interesting effect that does this is the spacetime shear. It is intriguing to find that the physical agencies such 
as the spacetime shear, and related inhomogeneities in matter density distribution within a dynamically collapsing 
cloud, could naturally delay the formation of trapped surfaces during gravitational collapse (39j . (40| . In other words, 
such physical factors do naturally give rise to naked singularity phases in the collapse where the formation of apparent 
horizon and the trapped surfaces is delayed. Even though the matter densities are arbitrarily large and growing, the 
shear could distort the trapped surface geometry in such a manner so as to avoid the trapping of light and facilitates 
the escape of null rays from such ultra-dense regions. 
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